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Abstract. The search for MSSM Higgs bosons will be an important goal at the LHC. In order to 
analyze the search reach of the CMS experiment for the heavy neutral MSSM Higgs bosons, we 
combine the latest results for the CMS experimental sensitivities based on full simulation studies 
with state-of-the-art theoretical predictions of MSSM Higgs-boson properties. The experimental 
analyses are done assuming an integrated luminosity of 30 or 60 fb~^. The results are interpreted 
as 5 a discovery contours in MSSM Ma -tan l3 benchmark scenarios. Special emphasis is put on the 
variation of the Higgs mixing parameter ^. While the variation of /i can shift the prospective dis- 
covery reach (and correspondingly the "LHC wedge" region) by about Zitan/3 = 10, the discovery 
reach is rather stable with respect to the impact of other supersymmetric parameters. Within the 
discovery region we analyze the accuracy with which the masses of the heavy neutral Higgs bosons 
can be determined. An accuracy of 1-4% should be achievable, depending on Ma and tan/3. 

PACS. 14.80.Cp Non-standard-model Higgs bosons - 12.60.Jv Supersymmetric models 



1 Introduction 

Identifying the mechanism of electroweak symmetry 
breaking will be one of the main goals of the LHC. The 
most popular models are the Higgs mechanism within 
the Standard Model (SM) and within the Minimal Su- 
persymmetric Standard Model (MSSM) [Ij. Contrary 
to the case of the SM, in the MSSM two Higgs doublets 
are required. This results in five physical Higgs bosons 
instead of the single Higgs boson of the SM. These are 
the light and heavy CP-even Higgs bosons, h and H, 
the CP-odd Higgs boson. A, and the charged Higgs bo- 
son, H^. The Higgs sector of the MSSM can be spec- 
ified at lowest order in terms of the gauge couplings, 
the ratio of the two Higgs vacuum expectation val- 
ues, tan/3 = V2/V1, and the mass of the CP-odd Higgs 
boson. Ma- Consequently, the masses of the CP-even 
neutral Higgs bosons and the charged Higgs boson are 
dependent quantities that can be predicted in terms of 
the Higgs-sector parameters. Higgs-phenomenology in 
the MSSM is strongly affected by higher-order correc- 
tions, in particular from the sector of the third gener- 
ation quarks and squarks, so that the dependencies on 
various other MSSM parameters can be important. 

The current exclusion bounds within the MSSM [21 
[31U and the prospective sensitivities at the LHC are 
usually displayed in terms of the parameters Ma and 
tan /3 that characterize the MSSM Higgs sector at low- 
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est order. The other MSSM parameters are convention- 
ally fixed according to certain benchmark scenarios O 
ijj . We focus here [7J on the 5 a discovery contours for 
heavy MSSM Higgs bosons, i.e. the lower bound of the 
"LHC wedge" , within the "m™'^'^ scenario" . For the in- 
terpretation of the exclusion bounds and prospective 
discovery contours in the benchmark scenarios it is im- 
portant to assess how sensitively the results depend on 
those parameters that have been fixed according to the 
benchmark prescriptions. Consequently, we investigate 
how the 5 a discovery regions in the M^^-tan (3 plane 
for the heavy neutral MSSM Higgs bosons obtainable 
with the CMS experiment at the LHC depend on the 
other MSSM parameters. 



2 The analysis 

The search for the heavy neutral MSSM Higgs bosons 
at the LHC will mainly be pursued in the b quark as- 
sociated production with a subsequent decay to r lep- 
tons |8','9^,'T0] . In the region of large tan (3 this produc- 
tion process benefits from an enhancement factor of 
tan^ (3 compared to the SM case. The main search 
channels arq3 (here and in the following (p denotes the 

^ In our analysis we do not consider diffractive Higgs 
production, pp — > p(BH(Bp [11]. For a detailed discussion of 
the search reach for the heavy neutral MSSM Higgs bosons 
in diffractive Higgs production we refer to Ref. [12] . 
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two heavy neutral MSSM Higgs bosons, (f> = H,A): 

hh(t),4>^ T+T- 2 jets (1) 
bb(j), (f) T+T^ ^ ^ + jet (2) 
bb(j), (f) T^T^ ^ e + jet (3) 

The analyses were performed with full CMS de- 
tector simulation and reconstruction for the follow- 
ing three final states of di-r-lepton decays: t+t~ — » 
jets [13], r+T^ e -I- jet [14; and r+r^ ^ /i + 
jet [T3]. The Higgs-boson production in association 
with b quarks, pp —> bb(j), has been selected using single 
6-jet tagging in the experimental analysis. The kine- 
matics of the gg bbcj) production process (2 3) 
was generated with PYTHIA [16]. The backgrounds 
considered in the analysis were QCD muli-jet events 
(for the TT jets mode), ti, 55, Drell-Yan produc- 
tion of Z, 7*, ly-l-jet, Wt and Trbb. AU background 
processes were generated using PYTHIA, except for 
T~^T~bb, which was generated using CompHEP [T7] . 



4> T+T- jets, 60 fb"^ 


Ma [GeV] 


200 


500 


800 


Ns 


63 


35 


17 


Scxp 


2.5 X 10"* 


2.4 X 10"^ 


3.6 X 10"^ 




0.176 


0.171 


0.187 


AM^/M^ [%] 


2.2 


2.8 


4.5 



Table 1. Required number of signal events, Ns, with 
£ = 60 fb~^ for a 5 cr discovery in the channel — > t'^t~ 
jets. Furthermore given are the total experimental selec- 
tion efficiency, Ecxp, the ratio of the di-r mass resolution to 
the Higgs-boson mass, Rm^, and the expected precision of 
the Higgs-boson mass measurement, AM^/M^, obtainable 
from Ns signal events. 



(j) T+T- ^ e + jet, 30 fb"^ 


Ma [GeV] 


200 


300 


500 


Ns 


72.9 


45.5 


32.8 




3.0 X 10"^ 


6.4 X 10"^ 


1.0 X 10"^ 


Rm^ 


0.216 


0.214 


0.230 


AM^/M^ [%] 


2.5 


3.2 


4.0 



Table 2. Required number of signal events, Ns, with 
£ = 30 fb~^ for a 5 cr discovery in the channel (f> t'^t~ — + 
e 4- jet. The other quantities are defined as in Tab.[T] 



(j> ^ T+T' n+ jet, 30 fb~^ 


Ma [GeV] 


200 


500 


Ns 


79 


57 


£'cxp 


7.0 X 10~^ 


2.0 X 10"^ 


Rm^ 


0.210 


0.200 


AM^/M^ [%] 


2.4 


2.6 



Table 3. Required number of signal events, Ns, with 
£ = 30 fb~^ for a 5 cr discovery in the channel — > t+t~ 
H + jet. The other quantities are defined as in Tab. [1] 



The results quoted in Tabs. [Tj- [3] for the required 
number of signal events depend only on the Higgs- 
boson mass, i.e. the event kinematics, but are indepen- 
dent of any specific MSSM scenario. In order to deter- 
mine the 5 (7 discovery contours in the M^i-tan (3 plane 
these results have to be confronted with the MSSM 
predictions. The number of signal events, iVov, for a 
given parameter point is evaluated via 

TVcv = X cr^^^ X BR{(j} ^ r+r") X BR^^ x e^^p . (4) 

Here C denotes the luminosity collected with the CMS 
detector, aj^i^ is the Higgs-boson production cross sec- 
tion, BR(0 — ^ T+T') is the branching ratio of the 
Higgs boson to t leptons, BRt-,- is the product of the 
branching ratios of the two r leptons into their respec- 
tive final state, 

BR(r ^ jet + a:) w 0.65 , (5) 
BR(t ^ ^ + A) « BR(r ^ e + A) « 0.175 , (6) 

and £oxp denotes the total experimental selection effi- 
ciency for the respective process (as given in Tabs. [Tj- 
[3]) . For our numerical predictions of total cross sections 
(see Ref. [TS] and references therein) and branching ra- 
tions of the MSSM Higgs bosons we use the program 
FeynHiggs 19,20,21.22 . We take into account effects 
from higher-order corrections and from decays of the 
heavy Higgs bosons into supersymmetric particles. 

In spite of the escaping neutrinos, the Higgs-boson 
mass can be reconstructed in the H, A —> tt channel 
from the visible r momenta (r jets) and the missing 
transverse energy, E!^^^^, using the coUinearity approx- 
imation for neutrinos from highly boosted r's. In the 
investigated region of Ma and tan/? the two states A 
and H are nearly mass-degenerate. For most values 
of the other MSSM parameters the mass difference of 
A and H is much smaller than the achievable mass 
resolution, and the difference in reconstructing the A 
or the H will have no relevant effect on the achiev- 
able accuracy in the mass determination. The preci- 
sion AM^/M^ shown in Tabs. [T]~[3]is derived for the 
border of the parameter space in which a 5 cr discov- 
ery can be claimed, i.e. with A^5 observed Higgs events. 
The statistical accuracy of the mass measurement has 
been evaluated via 

AM^IM^=RmJVn~s ■ (7) 

A higher precision can be achieved if more than Ng 
events are observed. The corresponding estimate for 
the precision is obtained by replacing Ns in eq. ([7]) by 
the number of observed signal events, A'ov It should 
be noted that the prospective accuracy obtained from 
eq. d?]) does not take into account the uncertainties of 
the jet and missing energy scales. In the t+t~ — > 
jets mode these effects can lead to an additional 3% 
uncertainty in the mass measurement _^13^. 

3 Numerical results for the LHC wedge 

We have evaluated A'ov in the rri^^^ benchmark sce- 
nario [5l[6] as a function of AIa and tan /3. For fixed Ma 
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we have varied tan f3 such that Ncv = Ns (as given in 
Tabs. [T]- 131). This tan/3 value is then identified as the 
point on the 5 a discovery contour corresponding to the 
chosen value of Ma- In this way we have determined 
the 5 fj discovery contours for the m™'*'^ scenario for 
H = ±200, ±1000 GeV.0 

In Fig. [T] we show the 5cr discovery contours ob- 
tained from the process 660, cj) t^t^ for the final 
states T+T^ jets, t+t~ — > e + jet and r+r" — > 
/X + jet. The 5 a discovery contours are affected by a 
change in fi in two ways. Higher-order contributions, 
in particular the ones associated with Z\f, [24! , modify 
the Higgs-boson production cross sections and decay 
branching ratios. Furthermore the mass eigenvalues of 
the charginos and neutralinos vary with ^, possibly 
opening up the decay channels of the Higgs bosons to 
supersymmetric particles, which reduces the branching 
ratio to r leptons. 

As expected from the discussion of the Z\f, correc- 
tions in Refs. [6l[7], the variation of the 5cr discovery 
contours with fj, can be sizable. In the t~^t^ jets 
channle (top plot in Fig. [T]) a shift up to A tan /? = 12 
can be observed for Ma = 800 GeV. For low Ma val- 
ues (corresponding also to lower tan/3 values on the 
discovery contours) the variation stays below A tan /3 = 
3. In the no- mixing scenario the variation does not ex- 
ceed A tan /3 = 5. The 5 a discovery regions are largest 
for /i = —1000 GeV and pushed to highest tan /3 values 
for /i = +200 GeV. In the low Ma region our discovery 
contours are very similar to those obtained in Ref . [5] . 
In the high Ma region. Ma ^ 800 GeV, correspond- 
ing to larger values of tan /3 on the discovery contours, 
our improved evaluation of the 5 a discovery contours 
gives rise to a shift towards higher tan/3 values com- 
pared to Ref. [5] of about A tan (3 = 8 (mostly due to 
the up-to-date experimental input). 

The results for the channel r+r^ e + jet are 
shown in the middle plot of Fig. [T] The resulting shift 
in tan /3 reaches up to A tan/3 = 8 for Ma = 500 GeV. 
Finally in the bottom plot of Fig. [T] the results for the 
channel r+r^ — *■ /i + jet are depicted. The level of 
variation of the 5 a discovery contours is the same as 
for the e + jet final stated 

In Ref. [7] it has been shown that the effects visible 
in Fig. [1] arising from the variation of /j, are a mixture 
of two effects: the change in the bottom Yukawa cou- 
pling via Ah and the impact on the heavy Higgs decay 
channels of possible additional decays to charginos or 
neutralinos. The variation of other parameters enter- 
ing the radiative corrections is comparably small. 
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Fig. 1. Variation of the 5a discovery contours obtained 
in the m™^" scenario for different values of jj, from the 
channels bb(j),(j) — > r+r^ jets (top), e + jet (middle), 
fi + jet (bottom). 



A corresponding analysis in benchmark scenarios 
fulfilling cold dark matter constraints can be found in 
Ref. [23]. 

^ Since the results of the experimental simulation for this 
channel are available only for two Ma values, the interpo- 
lation is a straight line. This may result in a slightly larger 
uncertainty of the results compared to the other two chan- 
nels. 



4 Numerical results for the Higgs-boson 
mass precision 

The expected statistical precision of the heavy Higgs- 
boson masses is evaluated according to eq. ([7]). In Fig. [2] 
we show the expected precision for the mass measure- 
ment achievable from the channel 660, </> —>■ t~^t~ using 
the final state r+r~ jets. Within the 5(t discovery 
region we have indicated contour lines corresponding 
to different values of the expected precision, AM/M. 
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The results are shown in the m™^^ benchmark scenario 
for II = -200 GeV (upper plot) and fi = +200 GeV 
(lower plot). We find that experimental precisions of 
AM^/M^ of 1-4% are reachable within the discovery 
region. A better precision is reached for larger tan /3 
and smaller Ma as a consequence of the higher num- 
ber of signal events in this region. The other channels 
and other values of ^ discussed above yield qualita- 
tively similar results to those shown in Fig. [2j 



C 47 
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Fig. 2. The statistical precision of the Higgs-boson 
mass measurement achievable from the channel bb(f>, (f> 
T^T^ — > jets in the m™™ benchmark scenario for fi = 
-200 GeV (top) and n = -1-200 GeV (bottom) is shown 
together with the 5 a discovery contour. 



Acknowledgements 

We thank S. Gennai, A. Kalinowski, R. Kinnunen and 
S. Lchti for collaboration on the work presented here. 
Work supported in part by the European Commu- 
nity's Marie-Curie Research Training Network under 
contract MRTN-CT-2006-035505 'Tools and Precision 
Calculations for Physics Discoveries at Colliders'. 



References 

1. H. Nilles, Phys. Kept. 110 (1984) 1; H. Haber and 
G. Kane, Phys. Kept. 117 (1985) 75; R. Barbieri, Riv. 



6 



10 



11 



12 



13 

14 
15 

16 

17 



18 



19 



20 



21 



22 



23 



24 



Nuovo Cira. 11 (1988) 1. 
[LEP Higgs working group], Eur. Phys. ,J. 
(2006) 547 [arXiv:hep-ex/0602042]. 
V. Abazov et al. [DO Collaboration] , Phys. Rev. Lett. 
95 (2005) 151801 [arXiv:hep-ex/0504018]; Phys. Rev. 
Lett. 97 (2006) 121802 [arXiv:hep-ex/0605009]; DO 
Note 5331-CONF. 

A. Abulencia et al. [CDF Collaboration], Phys. Rev. 
Lett. 96 (2006) 011802 [arXiv:hep-ex/0508051]; CDF 
note 8676. 

M. Carena, S. Heinemeyer, C. Wagner and G. Wei- 
glein, Eur. Phys. J. C 26 (2003) 601 [arXivihep- 
ph/0202167]. 

M. Carena, S. Heinemeyer, C. Wagner and G. Wei- 
glein, Eur. Phys. J. C 45 (2006) 797 [arXivihep- 
ph/05 11023]. 

7. S. Gennai, S. Heinemeyer, A. Kalinowski, R. Kin- 
nunen, S. Lehti, A. Nikitenko and G. Weiglein, to 
appear in Eur. Phys. J. C, arXiv:0704.0619 [hep-ph]. 

8. ATLAS Collaboration, Detector and Physics Perfor- 
mance Technical Design Report, CERN/LHCC/99-15 
(1999). 

9. K. Cranmer, Y. Fang, B. Mellado, S. Paganis, 
W. Quayle and S. Wu, hep-ph/0401148. 
CMS Physics Technical Design Report, Volume 2. 
CERN/LHCC 2006-021, see: 
cmsdoc.cern.ch/cms/cpt/tdr/ . 
M. Albrow and A. Rostovtsev, [arXiv:hep- 
ph/0009336]; V. Khoze, A. Martin and M. Ryskin, 
Eur. Phys. J. C 23 (2002) 311 [arXiv:hep- 
ph/0111078]; A. De Roeck, V. Khoze, A. Martin, 
R. Orava and M. Ryskin, Eur. Phys. J. C 25 
(2002) 391 [arXiv:hep-ph/0207042]; B. Cox, AIP 
Conf. Proc. 753 (2005) 103, [arXiv:hep-ph/0409144]; 
J. Forshaw, arXiv;hep-ph/0508274. 
S. Heinemeyer, V. Khoze, M. Ryskin, W. Stirling, 
M. Tasevsky and G. Weiglein, arXiv:0708.3052 [hep- 
ph]. 

S. Gennai, A. Nikitenko and L. Wendland, CMS Note 
2006/126. 

R. Kinnunen and S. Lehti, CMS Note 2006/075. 
A. Kalinowski, M. Konecki and D. Kotlinski, CMS 
Note 2006/105. 

T. Sjostrand et al., Comput. Phys. Commun. 135 
(2001) 238 [arXiv:hep-ph/0010017]. 
E. Boos et al. [CompHEP Collaboration], Nucl. 
Instrum. Meth. A 534 (2004) 250 [arXiv:hep- 
ph/0403113]. 

T. Hahn, S. Heinemeyer, F. Maltoni, G. Weiglein and 
S. Willenbrock, arXiv:hep-ph/0607308. 
S. Heinemeyer, W. HoUik and G. Weiglein, Com- 
put. Phys. Commun. 124 (2000) 76, [arXiv:hep- 
ph/9812320]; see: www.feynhiggs.de . 
S. Heinemeyer, W. Hollik and G. Weiglein, Eur. Phys. 
J. C 9 (1999) 343 [arXiv:hep-ph/9812472]. 

G. Degrassi, S. Heinemeyer, W. Hollik, P. Slavich 
and G. Weiglein, Eur. Phys. J. C 28 (2003) 133 
[arXiv:hep-ph/0212020]. 

M. Frank, T. Hahn, S. Heinemeyer, W. Hollik, 

H. Rzehak and G. Weiglein, JHEP 0702 (2007) 047 
[arXiv : hep-ph /0611326]. 

J. Ellis, T. Hahn, S. Heinemeyer, K. Ohve and G. Wei- 
glein, arXiv:0709.0098 [hep-ph]. 

M. Carena, D. Garcia, U. Nierste and C. Wag- 
ner, Nucl. Phys. B 577 (2000) 577 [arXiv:hep- 
ph/9912516]. 



